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RADIATION FROM INTENSE ELECTRON BEAMS ASSOCIATED WITH 

THE CERENKOV MECHANISM 
John R. Neighbours and Fred R. Buskirk 

Abstract 

This report summarizes the Cerenkov radiation to be expected 
from an intense electron beam for both the coherent and incoherent 
modes of radiation . The effects of a finite length beam are dis- 
cussed with regard to oeam diffraction, radiated power, and bunch 
snape. The appendices contain a discussion of the beam radiation 
pattern. 
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RADIATION FROM INTENSE ELECTRON BEAMS ASSOCIATED WITH 



THE CERENKOV MECHANISM 

1 . INTRODUCTION 

Normally a charge or system of charges moving at constant 
velocity does not radiate electromagnetic energy, but if the 
velocity of the charge is greater than the velocity of light in a 
medium, radiation occurs. It is broad band; the intensity 
increases with freauency; and the propagation is confined to a 
cone, much like a shock wave in air. Equations (?) and (3) in a 
later section describe these results. In previous work 1 “ 4 , we 
described theory and experiments showing that this mechanism 
produces microwaves from bunched electron beams, in addition to 
the usual optical radiation, as studied by Cerenkov and many 
others . 

In this report we consider aspects of coherent Cerenkov 
radiation, with emphasis on radio and microwave frequencies. As 
will be shown, power levels may become very high, in the megawatt 
range, so that serious consideration should be given to this 
mechanism. The scope of this report is to outline the systematics 
of a theory in which the dielectric (v/hich lowers the velocity of 
light) is considered to be linear and non-dispersive . The aspects 
considered which are not treated in the extensive literature are 
calculations for electrons occurring in finite size bunches, both 
single and repeated. Future work should aim at the effect of 
ionization in the dielectric medium, betatron oscillations of the 
beam electrons, and other mechanisms which could seriously modify 
or degrade the radiation. 
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SUMMARY OF PREVIOUS WORK 



In the earlier work 1 - 4 , it was assumed that electrons were 
emitted from an accelerator in bunches, described by a charge 
density p(r), and all had the same velocity v such that the bunch 
did not change shape in time. The current density is then given 
by 

J(?) = P v (1) 



A. Intensity 

Under these conditions, the Cerenkov radiation in a medium 
with the velocity of radiation c < v is 



d_E 

dx 



dE 

dx 



iod uj sin 2 0 q 2 F 2 (k) 
c 

o)o) 0 sin 2 0 c q 2 F 2 (k) 



(2) 

(3) 



Eq. (2) refers to a single bunch, the radiation has a 
continuous frequency spectrum, and the formula gives the energy 
radiated at angular frequency o) in a range do), per unit length. 
Eq. (3) is the energy-per-unit path per bunch for the case of 
periodic bunches, emitted at an angular frequency o) Q . Here o) is 
a harmonic of o)q, in contrast to the continuous frequency 
spectrum of Eq. (2). In both cases, q is the total charge of one 
bunch, and F(>c) is a form factor for the charge distribution 

F ( k ) = q 1 / / / p ( r ) e l] ^’ r dV (4) 

0 C is the Cerenkov angle, 
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q is the charge per bunch, and k is propagation vector for the 
emitted radiation, such that k = to/c, and the direction is along 
the Cerenkov cone. 

B. Polarization 

The electron field is in the plane of the observer and the 
electron path, as it is for usual optical Cerenkov radiation. 

C. Coherence 

If the charcre is a point charge, F = 1 and in fact, for a 
finite bunch, we will have F ~ 1 for k£ << 1, where k = 2n/\. 

This corresponds to X > Z, where Z is the size parameter of the 
bunch. For these long wave lengths, either formula, Eq. (2) or 
(3), shows that the intensity is found by setting q = ne where n 
is the number of electrons in the bunch. Thus the intensity is 
proportional to n 2 e 2 . As higher frequencies are examined, the 
factor F becomes small, and finally, for very high frequencies 
(optical range), we would expect incoherent radiation, in which q 2 
is replaced by ne 2 . 

D. Diffraction Effect 

An important result is that if the region of emission of 
radiation is finite, of length L, the emitted radiation is no 
longer confined to the Cerenkov cone 9 = Q c but is spread about 
that angle and the spreading depends on the length of the path in 
the gas and the wave length of the radiation. It is more closely 
related to the diffraction pattern of an end fire antenna, rather 
than the broadside array; the latter is related to the single 
slit of optics. This effect is independent of bunch size and 
repetition rate. Associated with diffraction is an enhancement, 
qualitatively described below. 
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E. Enhancement from Diffraction 



In most diffraction or interference situations in optics, 
the energy is usually redistributed in angle if parameters are 
changed. For example, if the width of a single slit is changed 
but the total power passing through the slit remains constant, the 
total power integrated over the diffraction pattern on a screen 
remains constant. In the Cerenkov case, the development of Ref. 3 
showed that for small Cerenkov angles, the intensity at 9 C 
remains constant as L decreases, but the total power radiated at 
angles 0 > 9 C increases dramatically. The total power radiated 
increases by factors of up to 50 in our experimental situation. 

F. Form Factor Effect 

The basic Cerenkov radiation strength tends to increase 
with frequency, for either coherent or incoherent emission. 
Starting at low freauencies, the coherent emission increases with 
w until F becomes small, at u = u^. At this frequency, the wave 
length is short enough so that radiation emitted from electrons at 
the front and rear of the bunch begin to interfere destructively. 
This leads to destructive interference and lowered emissions of 
power above . Experiments in principle could demonstrate this 
effect, but at this time we have been unable to make measurements 
in the high frequency microwave range ( X ~ 2 mm) which would 
quantitatively demonstrate the predictions . 

As a further remark, it should be possible with precise 
measurements to determine size parameters of the electron bunches 
at a distance, or, in other words, to do beam diagnostics by this 
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method. These possibilities will be explored in another report. 
3. Approximate Results for Radiated Power 

The results mentioned previously may be related to more 
conventional parameters of pulse length and peak current. If we 
assume a pulse in the form of a uniform cylinder, J = pv and the 
current I = JA, q = pA£ = pAvT where SL = the length of the bunch 

and T = the bunch length measured as a time. Thus we find 

q = IT (6) 

For a Gaussian pulse, the charge density is 

-(x 2 + y 2 )/a 2 -b 2 z 2 /b 2 (7) 

p = p 0 e e 

It may be shown that the peak current in the middle of the 
pulse (at z = 0) is 



I = 
o 



ifb 



and that the normalization constant 
charge of the bunch by 
q = P, 



3/2 2 

it a b 



( 8 ) 

0 Q is related to the total 

(9) 



Then we may let the effective length of the bunch be l = bu 

and use q = I 0 T = I 0 l/v , just as for the uniform cylinder. 

Now consider the sum over frequencies in Eq. (3). Assuming 

that F(k) is constant up to some k which corresponds to 

max r 

u = Nii) n , we find 
m u 



o 



£<dw 0 = £n a > 0 2 



= N ( N+l ) a) 0 2 /2 

“ N 2 u > 0 2 /2 

= <%i 2 /2 



Thus 





(10) 



For a single bunch, the sum over the continuous frequency 
distribution (2) gives 



Thus the results of Eqs . (10) and (11) above may be used for 

either single or repeated pulses to calculate the total power 
radiated at all frequencies. 

Now we may use the results of Eq . (6)- (9) to re-express the 
results of Eq. (?) or (3). The latter may be multiplied by v to 
convert to power, then summed or integrated over o> for the total 
power emitted, which is 



This holds for either periodic or a single bunch. In either 
case it is the rate at which electromagnetic radiant energy is 
emitted while the bunch is in flight, per bunch (in the case of 
periodic bunches). 

Noting that f m = l/T holds, which is equivalent to the 
uncertainty relation in quantum theory or the bandwidth-pulse 
length relation in communications, 




(ID 




2 2 2 

sin 2 0 I n T f 
c 0 m 



( 12 ) 
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p = JLiLa- v sin 2 0 I 2 
■» c 0 



( 13 ) 



Thus the power emitted for coherent radiation is proportional 
to the square of the beam current. This makes explicit one of the 
main points of this report: that the emitted power increases 
rapidly with the beam current. Thus, although weak beams emit 
little radiation, the radiation becomes important for large 
currents. When we note that the radiation emitted is proportional 
to Iq 2 but the energy in the beam is proportional to I Q , so that 
relative efficiency as a radiator increases in proportion to I Q . 
This may be made quantitative by multiplying (13) by L/v to get 
the total energy emitted by the bunch during its flight through a 
path L, where L is the range. The energy put into the bunch is 
IqT, so that the efficiency £, defined as the ratio of the energy 
radiated (Cerenkov) to the total energy of the beam becomes 



where V Q is the energy of the beam expressed in volts and L is the 
range of the beam in air. Strictly speakinq, L is the distance at 
which Cerenkov radiation ceases. The precise value of L would 
depend on the relative importance of bremsstrahlung, ionization, 
and instabilities which could disperse the electron bunch, and 
these aspects should be considered in greater detail . 

The overall systematics of the radiation associated with the 
Cerenkov mechanism are expressed by Eq. (12), (13), and (14). The 

following features should be noted: 



e 




(14) 
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A. The power radiated per bunch (eq. 12, 13) is proportional 

to the square of the peak beam current, I Q 2 , and is independent of 
the pulse length T. This depends on the approximation that 

f m T = 1, where T is the pulse length and f m is the maximum 
frequency for which the fourier transform of the pulse has 
significant strength. This expression for the power radiated is 
the energy-per-unit time while the bunch is in flight, from the 
accelerator to the point where it ceases to radiate effectively. 
The latter point in the simplest consideration is the point at 
which the bunch has lost enough energy so that it is below the 
threshold for Cerenkov radiation, which is 20 MeV for electrons in 
normal air . 

B. As T is decreased, assuming I Q is constant, the power 
radiated remains constant. But the frequency at which most of the 
radiation occurs is increased as T is decreased. 

C . The efficiency (Eq. 14) depends on I 0 /T showing that for 
high currents and short bunches, the Cerenkov mechanism could not 
only lead to sizable amounts of radiation but also dominate the 
other modes of energy loss. 

D . As a guide to what may be expected in the case of plasma 
shielding, we might expect that some length l 1 of the nose of the 
beam bunch might be unshielded by the plasma. From the above 
considerations, the length i ' would not change the total power 
radiated, but would change dramatically the freauency at which the 
radiation occurs. This is because if only the nose region of 
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of length l' radiated, the total effective radiating charge would 
be proportional to But the frequency of the maximum radiation 

would increase as 1/2.', so the power radiated would remain 
constant as l' is varied. 
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APPENDIX A - EXTREME VALUES OF THE RADIATION FUNCTION 



The radiated power per unit solid angle W(v,n) can be written as the 
product of two factors: 

W = Q R 2 (Al) 

where Q is a function having the dimensions of power (watts) and R is a 

dimensionless radiation function, 
ye V 0 2 

Q = q 2 ( A2) 

8n 2 

R = (kL) SinO F(k) I(u) (A3) 

Here, y is the permeability of the medium, c is the velocity of light in 
the medium, and q is the charge of an individual electron bunch. If an RF 
Linac is the source of the electrons, the frequency v Q is that of the Linac ; 
however, in all cases, v Q corresponds to the spatial periodicity, l - c/v Q , 
of the electron bunches. 

In (3), k is the wave vector of the emitted radiation, L is the beam 
interaction length, and F(k) is the form factor defined such that 



£ 0 ’(k) = q F(k) 

and £Q*(k) is the Fourier transform of the charge distribution of an 

charge bunch. The diffraction function 
j _ Sin u 
u 

is a function of the diffraction variable 

VT \ 

u= (Cos0 -CosOl 
~ C 



( A4) 



individual 



(A5) 



(A6) 



and consequently, R is a somewhat comnlicated function of 0 and kL . 
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MINIMUM VALUES 



By inspection, it is obvious that and thus W will have minima when 
l(u) = 0. These true zeroes of the function occur when u is equal to integral 
values of tt; i.e., u = m tt. Substituting this condition in (A6) leads to the 

angular values for the zeroes in the radiated power 

Cos0=CosO - m -p— (A7) 

which can be written in terms of the wavelength (in the medium) of the 
radiation as 

Cos0=Cos0-mf^} (A8) 

This relation is exact and gives the zeroes of R for all allowed values of 
m for which Cos 0 < 1. I(u) is a maximum at m = 0 (which gives u = 0), and the 

maximum in R tends to occur somewhere in that vicinity. The first zero is that 

one which is adjacent to the maximum and is given by (A8) when m = +1. Values 
of this angle as a function of are shown as a Log -Log plot in Fig. 1. The 
relationship between 0 and is approximately that of the inverse square root 
over the range 10 < L/ X < 100. Within this range, the angular value (in 
degrees) of the first minimum is given by the empirical relation 

01 = - ? 8 - 2 - - ■■ ( A.9 ) 

VL/X 

The physical range of 0 is from 0 to tt. Substitution of the successively 
larger positive values of m = 2, 3, 4... into (A7) or (A8) gives a set of 
increasing angles for the functional zeroes. If X/L « 1, many zeroes will 
occur before the physical limit is reached. 

In addition, if X/L << 1, there is the possibility of having functional 
zeroes at angular values less than that of the main peak. These occur with the 
substitution of successive negative m values into (A7) or (A8) , giving rise to 
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subsidiary radiation maxima between the beam line and the main diffraction 
peak. 

MAXIMUM VALUES 

To find the maximum values of ^ one proceeds in the usual fashion by 
taking the derivative and setting it to zero. 



aw = 2 o R = 0 

Jo 2 Q R Jo 0 



(A10) 



3R 



This equation is satisfied for either R = 0 or -r-r = 0. The first condition is 

o d 

that considered above for the zeroes of the function. The remaining one, 

3R 

-g-g = 0, pertains to the maxima. Thus, it is sufficient to consider maximum 
values of R which will also be maximum values for W ~ R^. 

From (A3), the derivative of R is 

H - kL { sin9 [f + I 4S] + FI CosQ } = 0 (All) 



and the derivative of the diffraction function I(u) can be written 
31 dl 3u Tt kL . 

VS dt 15 1 — Sln0 



(A12) 



Substituting this value in (All) and discounting the solution for kL = 0 gives 



the condition for maxima 



{ SinO [i' F Sin9 + I -||] + FI CosO} = 0 



(A13) 



In order to proceed further, the functional dependence of the form factor 
F must be known. 
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A gausian form factor has been assumed in this and other reports. 
F(k) = Exp (-a 2 Cos 2 0) 

where a= kb/4 corresponding to a one-dimensional charge distribution 

Pq' (r) = q Exp (- — ) 

b 2 



( A 1 4 ) 



(A15) 



The derivative of this form factor is 
3F 9 

10 = a Sin 29 F 



(A16) 



which is small for the parameters of interest in Linac experiments; i.e., 
b = 0.24 cm, microwave frequencies, and angles up to approximately 25°. 
Therefore, a first approximation to solving (A13) can be obtained by setting 
the derivative of the form factor in (A13) to zero. The resulting aDproximate 
maxima condition is independent of F and can be rearranged so that I and its 



derivative are on the same side. 

- X< = _L_ .Cos 9 
I(u) kL Sin 2 9 

The ratio of I 1 to I can be expanded as a power series in u as 



V_ 

I 



u u 
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(A17) 



( A18) 



3 43 

which converges rapidly for small values of u. Substituting the first term in 
the expansion into (A17) and subsequently eliminating the diffraction variable 
through its definition (A6) leads to an approximate equation for the maxima, 
which is cubic in y = Cos9. 



-Y 3 ~ Cos 3 Y 2 “ [l + —1 Y + Cos 9 = 0 

^ 2 C 



n 

x ‘ 



( A19) 



In this equation, the beam parameters enter through the charge velocity in 
Cos 9 C = c/v and the radiation wavelengths in x = kL = 2 tt L/ A. Numerical 
solutions can be carried out on a hand calculator with results within 
approximately 0.3° of those obtained from large-scale computer calculations. 
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Results for the major radiation lobe (the first maximum) are also plotted in 
Fig. 1. These also show a straight line dependence on a Log-Log plot and 
therefore can also be represented as a power law 

0 

m 

which is accurate to approximately 5%. 
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APPENDIX B - EFFECT OF BEAM LENGTH ON THE RADIATION FUNCTION 



In previous reports, several calculated curves of a radiation function 
D( 0 ) differs °nly bv constants from R 2 were presented. Fig. 2 is 

reproduced from Fig. 5 of Reference 4 and relabeled to be consistent with the 
notation of Appendix A. Fig. 2 shows the vartiation of the radiation function 
for 8.55 GHz radiation from the NPS Linac for beam interaction lengths in the 
vicinity of 100 cm. As noted previously in this and other reprots, increasing 
L tends to move the main diffraction lobe toward the beam line, and Appendix A 
gives an approximate formula for calculating the position of the first lobe. 
Fig. 3 (to different scales) shows this effect for larger beam paths from an S- 
Band Linac. As L increases, the principle lobe moves in toward the Cerenkov 
angle 9 C ; i.e., + 0 C , and the radiated power per unit solid angle at 

the maximum increases greatly. 



i 
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FIGURE CAPTIONS 



Fig;. 1 Cerenkov radiation extremum angles in degrees vs L/X. The wavelength X 
of the Cerenkov radiation is measured in the medium in which the path 
length is L. Curve B is for the first minimum, the solution to (A8) 
with m = +1. Curve A is for the first maximum, the solution to the 
approximate equation (A19). Over the range shown, A is straight and B 
is slightly curved. Equations for the curves are given in the text. 

Fig. 2 Calculated third harmonic (8.55 GHz) radiation patterns produced by 

electron bunches from an S-band (2.85 GHz) Linac . The path lengths in 
air are 70, 90, and 150 cm. Vertical scale is arbitrary and the same 
for all three curves. 

Fig. 3 Calculated third harmonic radiation patterns produced by electron 

bunches from an S-band Linac. The solid curve is for a path in air of 
1000 cm, and the dotted curve is for a path length of 10,000 cm. 
Vertical scale is arbitrary and the same for both curves but much 
reduced from that of Fig. 2. 
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